Universality of the Heisenberg limit for estimates of random phase shifts 
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The Heisenberg hmit traditionahy provides a lower bound on the phase uncertainty scaling as 
1/{N), where (N) is the mean number of photons in the probe. However, this limit has a number of 
loopholes which potentially might be exploited, to achieve measurements with even greater accuracy 
Here we close these loopholes by proving a completely rigorous form of the Heisenberg limit for the 
average error over all phase shifts, applicable to any estimate of a completely unknown phase shift. 
Our result gives the first completely general, constraint-free and non-asymptotic statement of the 
Heisenberg limit. It holds for all phase estimation schemes, including multiple passes, nonlinear 
phase shifts, multimode probes, and arbitrary measurements. 
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Accurate estimation of optical phase is a highly im- 
portant goal in quantum metrology. For example, many 
optical measurement schemes, ranging from temperature 
sensing to gravitational wave detection, rely on deter- 
mining the phase shift induced by an environmental pa- 
rameter [1 . Further, phase-based optical communication 
relies on being able to decode information carried in the 
phase of an optical field [2 . It is therefore of fundamen- 
tal theoretical and experimental interest to determine the 
maximum phase sensitivity that is possible when a given 
resource, such as the average photon number of a probe 
state, is available [3UTQ]. 

A key concept in quantum phase estimation is the 
Heisenberg limit [8^, TT, which is generally understood 
as the asymptotic scaling relation 

a{^)>k/{N), (1) 

for any proportional measure of phase resolution, cr(4>), 
where 4> is the estimate of the true phase ^. Here k is 
a constant of order one, and (N) is the average of the 
number operator which generates the phase shift, via the 
unitary operator ex.p{—iN^). In the simple case of a lin- 
ear, single-pass phase shift, N is the total photon number 
(across all modes which experience the phase shift). 

It is well known that the Heisenberg limit is valid 
for canonical phase measurements on single-mode probe 
fields [12j, and hence for covariant phase estimates on 
such fields However, more generally, the asymp- 

totic bound in Eq. ([T]) is merely heuristic and open to 
challenge. For example, general schemes for estimating 
phase shift parameters may include the use of multimode 
fields, multiple passes of probe states, nonlinear phase 
shifts, noncovariant phase measurements, and/or entan- 
gling joint measurements. 

Thus, there remain many loopholes which might poten- 
tially be exploited to achieve far greater phase measure- 
ment accuracies - for example, with asymptotic scalings 
better than 1/{N) and/or arbitrarily small values of k 
p^lHlGj . It is, therefore, of considerable importance to 



determine whether this is possible. 

In this paper all such loopholes are closed by giving 
a completely general constraint-free form of the Heisen- 
berg limit, valid for all phase estimation schemes in- 
cluding multiple passes, nonlinear phase shifts, multi- 
mode probes, and arbitrary measurements. Further, they 
are non-asymptotic, providing rigorous lower bounds on 
phase resolution valid for all values of (N). 

It will be shown analytically that 

6^>kA/{N+l), (2) 

where kA '= \/27r/e^ ~ 0.559 and (5^)'^ denotes the av- 
erage, over all applied phase shifts, of the mean square 
deviation of the phase estimate from the actual phase 
shift. The asymptotic form of this bound immediately 
yields the Heisenberg limit in Eq. ([T]), with k = kA- 
Moreover, numerical work strongly supports the conjec- 
ture that kA in Eq. ([2| can be replaced by kc ^ 1.376, 
which is asymptotically optimal as we show. 

It immediately follows from Eq. Q that no scheme 
can do better, on average, than the Heisenberg limit. It is 
crucial to the derivation of this bound that the phase is 
a priori completely unknown, or random. That this is no 
impediment to entanglement-enhanced phase estimation 
was recently shown experimentally [9 . Without such a 
condition, there is no universally valid Heisenberg limit 
of the form of Eq. ([2|. In fact, there are some singular 
range restrictions for which perfect phase resolution is 
possible in principle, for a finite average photon number, 
as discussed further below. 

All of the results in this paper are proven using only the 
assumption that the eigenvalues of the generator N are 
nonnegative integers. In the optics context, this allows 
for the very general case of multimode (possibly entan- 
gled) states, with a phase shift generated by an optical 
nonlinearity [SI [8], and where the kth mode, with photon 
number operator Nk^ passes through the phase shift pk 
times [4, 7 , i.e., N = "^ZkPki^kY ^ where q and pk are 
integers. Note that, in the case of nonlinear or multipass 
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phase shifts, bounds scahng inversely with (A^) do not 
preclude a better scaling (e.g. inverse square) in terms 
of mean photon number, as many have noted [3 [8]. 

Our results apply to any measurement whatsoever, in- 
cluding adaptive measurements where each mode k is 
measured separately, but the result of measuring mode 
k is used to alter the dynamics of mode k^ for k' > k 
as in Refs. [T, ^ , either by controlling pj^i or by apply- 
ing (possibly nonlinear) auxiliary phase shifts on those 
modes. This can be shown by extending and specializing 
the methods of Ref. [TT] . 

Giovanetti et al. have very recently obtained a related 
result, valid for estimation of any shift parameter, which 
for the case of phase estimation yields a bound on the 
average of the root mean square deviations for any two 
distinct applied phase shifts (j) and (j)': ^(A0<l> + A^/^) > 
t^{\^^>)/{N), where := |0 - 0'|/(A^4> + A^/4>) - 1. 
The function Hi has a maximum of ~ 0.074 at A^^/ ~ 4.7. 
This yields a Heisenberg limit, in an averaged sense, 
with scaling constant k ^ 0.074 [18 , valid when the 
strong constraint A^^/ ^ 4.7 is satisfied or (N) is suffi- 
ciently large. In contrast, the bound in Eq. Q is entirely 
constraint-free and far stronger. 

To proceed, consider a general phase estimation sce- 
nario in which a probe state, described by some density 
operator po on an arbitrary Hilbert space, undergoes a 
phase shift (j) (the value of the random variable ^) to be- 
come p(f) := e~'^^^ poe^^^ . Here N is any operator with 
nonnegative integer eigenvalues. A detection method M, 
described by some probability operator measure (POM) 
{M^}, is then used to estimate (j). That is, the proba- 
bility density of estimating a value ^, for the case of an 
actual phase shift 0, is given hy p{(j)\(j)) = tr[M^ p^]. 

The overall performance of the estimate, for a priori 
unknown phase shifts, may be characterised by the con- 
centration about zero of the average probability density 

m = ^ fj d4>p{e+m. (3) 

of the random variable 6 := <l> — <l> G [— tt, tt] (see Fig. 1). 
Note that for a measure of concentration linear in p{6)^ 
such as (O^), this concentration is just the average of 
the individual concentrations of p(^|^) about (j). More 
generally, for a measure of concentration that is convex 
(concave), its value for p{0) is a lower (upper) bound on 
the average concentration. 

We require three known results. First, the average 
phase distribution, p{0) in Eq. ([3| can be formally gener- 
ated by a covariant phase estimate M for the same probe 
state [19] (with POM {e-^^^Moe*^^}, where Mq := 
^/^^''^e^^^M^e-^^^). Second, the phase distribu- 
tion for a covariant phase estimate can be obtained via a 
canonical phase estimate, on a state with the same num- 
ber distribution of N [131 120]. Third, the canonical phase 




FIG. 1: (Color online) For a good estimate of a given phase 
shift 0, the corresponding probability density p{(j)\(j)) must be 
well concentrated about the value (j) — (f). In the scenario 
depicted in (a), the estimate is better for cj) — than for 
the cases = 02 and (j) = (jy^ (in the latter case the distribu- 
tion is offset from the actual phase shift 03 ). The degrees of 
concentration can be compared more directly by translating 
them to the origin, as in (b). Thus, equivalently, a good es- 
timate of (j) corresponds to the distribution of ^ = — 0, i.e., 
Pct){0) = p(^+(/)|(/)), being well concentrated about 0. It follows 
that how well the estimate performs overall, for randomly ap- 
plied phase shifts, is characterised by the degree to which the 
average distribution, p{0) in Eq. (jsj, is concentrated about 
^ = (dashed curve). 

and number distributions of any probe state are identical 
to those of a corresponding single mode field [71 [21] . 

As a consequence of these three results, the average 
phase distribution for any phase estimate is formally 
identical to the canonical phase distribution of some sin- 
gle mode field with the same number distribution. This 
yields the Lemma: 

Lemma. Any bound on the concentration of the canoni- 
cal phase distribution of a single mode field, under some 
constraint C on the photon number distribution, is also a 
bound on the concentration of the average phase distri- 
bution p{0) of an arbitrary phase estimate for any shift 
generator N having nonnegative integer eigenvalues, pro- 
viding that the probe state satisfies the same constraint C 
with respect to the distribution of N. 

Note that N may be any generator with nonnegative 
integer eigenvalues, and its spectrum may be a strict sub- 
set of the nonnegative integers (as in the case of optical 
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nonlinearities) . It is clear that all bounds will still hold 
under such a constraint on the spectrum, though it means 
that there may be tighter bounds. 

We now use the above Lemma to obtain the bound in 
Eq. ([2|. The average of the mean square deviation of an 
estimated phase shift from the actual phase shift, over 
all possible phase shifts ^, can be written as 



1.75 



(J4>)2 := (6") 



dO 



(4) 

where Var0<l> := J^^^ d^{^ - (j))^ p{^\^) [13 . The en- 
tropic uncertainty relation for the canonical phase and 
photon number observables of single mode fields ^13l [22] 
immediately yields, via the Lemma, H{Q) + H{N) > 
logg 27r for the entropies of O and TV. Further, for- 
mally extending p{6) from [— 7r,7r] to (—00,00), we ob- 
tain H{e) < I logj27re(^4>)2] [23j (where the inequality 
is strict since it is only saturated by Gaussian distribu- 
tions over the whole real line). Because entropy is max- 
imised for thermal distributions, H{N) < \og^{N + 1) + 
{N) logg(l + 1/(A^)). Noting the last term is always less 
than unity, these inequalities combine to yield 



27r 



(5) 



and, hence, Eq. ^ as desired. 

Thus we have shown that the Heisenberg limit holds 
in general. There is no way to take advantage of special 
measurement schemes in order to obtain measurements 
with scaling better than 1/(A^), provided the phase is 
initially unknown. By using (A^ + 1) in the denominator, 
rather than (A/"), we provide a bound that holds for all 
(N); it is clear that there cannot be bound of the form 
k/ (N) which is valid for all (TV), because the phase error 
cannot be larger than tt. Note that for the special case of 
m identical probe states, each with mean number (A^i), 
our lower bound evaluates to /c^/ {mNi + 1). 

A secondary question is whether, and by how much, 
the scaling constant in Eq. (|2| can be improved. Here we 
report strong numerical evidence for the conjecture that 



5^ > kc/{N^l), 



(6) 



with kc := 2{-za/S)^^^ ~ 1.37608, is the best possi- 
ble lower bound of this form, where za denotes the first 
(negative) zero of the Airy function. In particular, the 
minimum possible value of 6^ for any canonical phase es- 
timate on a single mode field, and therefore for any gen- 
eral phase estimate via the Lemma, may be numerically 
calculated as a function of (N) via standard techniques. 

In particular, one can use the method of Lagrange mul- 
tipliers, as described in Ref. [24j. In this case, rather than 
maximising (cos6>), one wishes to minimise (^^). Because 
{0'^) is symmetric and considered over the region [— 7r,7r], 
it can be expanded as a cosine series. Then the Lagrange 




FIG. 2: (Color online) The minimum possible value of 
{N + 1)^^, plotted as a function of {N) (upper curve). 
The aymptotic value is numerically indistinguishable from 
kc ^ 1.37608 in Eq. ^ (lower dotted curve). As per the 
text, the asymptotic value cannot be any less than kc- The 
intermediate curve replaces 6^ in Eq. ^ by the lower bound 
I — I cos^ + I cos 2^, for comparison purposes. This can be 
computed efficiently for larger values of (A), and shows the 
same asymptotic behaviour. 



multiplier problem may be formulated as an eigenvalue 
problem, and solved numerically. The resulting mini- 
mum values of are plotted as a product with (A' + 1) 
in the upper curve in Fig. 2. It is seen that the value of 
the product asymptotically approaches a minimum value 
from above that is numerically indistinguisable from kc- 



Moreover, because 



: cos -\- I cos 20 is no greater 



than the minimum expectation value of the first func- 
tion also places a lower bound on the minimum phase 
variance. This optimisation is computationally easier, 
and can be performed for larger values of (A"), because 
the eigenvalue problem obtained is sparse. The minimum 
values obtained are shown as the intermediate curve in 
Fig. 2. Again these results do not drop below /cc, pro- 
viding further evidence for the inequality ([6|. 

The constant kc can be obtained from the properties of 



the Holevo variance, {Sh^)^ 



1 [lailS]. For 



a canonical phase measurement on a single mode field, 
this has a tight asymptotic bound Sh^ ^ kc/{N)^ valid 
to first order in {N)~'^ [24l [261. The inequality chain 
|(e^®)p > (cos 6)2 > (1 - eV2)2 ^^^^ yig^g 



1= ^(m + 1) 



m>l 



for general estimates, via the Lemma. In the asymptotic 
limit only the m = 1 term contributes, to lowest order, 
and hence > kc/{N). Thus, the constant kc in the 
conjectured bound of Eq. (|6| is asymptotically optimal. 

Bounds on other direct measures of overall phase res- 
olution can be similarly obtained. For example, the 'en- 
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semble length' L{^) := e^^®\ which quantifies the effec- 
tive support length of p{0) [27^, has the analytic lower 
bounds L(l>) > 27re-^W > (27r/e)/(7V + 1), analogous 
to Eq. ([5|, using the same entropic uncertainty relation 
as before. The bound is stronger than, and implies, the 
bound in Eq. ([2|. The Lemma and general methods 
above can also be used to obtain similar bounds on over- 
all phase resolution for other constraints on A^, such as 
having a fixed maximum value nmax [28] . 

The above results establish the validity of the Heisen- 
berg limit, as a strong constraint on the overall per- 
formance of general phase detection schemes. In- 
deed, Eq. Q establishes a non-asymptotic analytic and 
constraint-free lower bound for average phase resolution. 
Furthermore, strong numerical support has been given 
for the conjecture that the numerator in Eq. Q can be 
replaced by the asymptotically optimal scaling constant 
/cc, as per Eq. ([6|. All results immediately generalise 
to bounds on the time resolution of any periodic system, 
with Hamiltonian H = Eq -\- hwN ^ via the replacement 
<i> uT. They can also be generalised to shift generators 
having arbitrarily negative integer eigenvalues [28l [29j . 

It follows that it is only possible to go beyond the 
Heisenberg limit when there is some restriction on the 
applied phase shift ^, i.e., when prior information about 
is available. In such cases one should not average 
over (j) with equal weighting, but instead should use a 
weighting corresponding to the prior probability. In- 
deed, there is the well known case where the phase shift 
is limited to values in the finite set {(j)k := 27rk/K}^ 
where k = 0,l,...i^ — 1 for some integer K > 1. 
Choosing the single mode probe state po = |?/^)(?/^|, 
with Itp) := K~^^'^ \n), and a phase estimate 

corresponding to measurement of the Hermitian opera- 



then the phase- 



tor 

shifted states are distinct eigenstates of ^k- Hence, 
the phase shifts can be resolved without error, even 
though the average photon number is finite. 

Another case is that where the error in the phase esti- 
mate is smaller than the Heisenberg limit only for some 
specific value of the phase. It might be imagined that 
one could take advantage of this measurement (without 
prior information about the phase) by using initial mea- 
surements to determine an initial estimate of the phase. 
Our results show that no such scheme can do better than 
the Heisenberg limit overall: when one takes into account 
the resources needed to initially estimate the phase, the 
complete measurement cannot beat the Heisenberg limit. 

It would be of interest to consider the case where such 
prior information is characterised by the properties of 
some prior probability density, ^(0), so that the factor 
l/(27r) in Eq. ([3| is replaced by integration over q{(l)). 
This case will be addressed in future work. 
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